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ABSTRACT 


The present work concerns itself with development of a software shell 
for fault diagnosis in rotating machinery installations. Expert 
system methodology of creating a knowledge-base and using inferencing 
techniques similar to that of a human expert are employed in 
developing the shell. Mathematical tools of Baye * s Conditional 
Probability Analysis and Influence Coefficient Matrix are used for 
the inferencing. 

The shell has been developed in GCLI5P. It is also referred to as 
knowledge-based system to indicate the fact it uses a stored 
knowledge-base (in the form of matrix of conditional probabilities) 
to decide on the order of the likely causes of vibration. The 
software lists the probable defects of the machinery in their order 
of likelihood of presence in the instal lat ion. The use of the shell 
has also been illustrated for the case of a High Speed Turbomachine 
in a Petroleum Plant. 
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CHAPTER 1 

INTRODUCTION 

Supervision of machinery health is a subject af increased importance 
in the course of progressive automation. Machinery condition 
monitoring is essential in order to anticipate problems in time in 
order to prevent complete failure. It serves the important task of 
ensuring operation within design constraints. Secondly, condition 
monitoring warns of changing or deteriorating conditions to prevent 
equipment failure and costly secondary damage. Early warning about 
the presence of a problem can be obtained while the problem is still 
localized. Thirdly condition monitoring provides an accurate 
description of machine condition in order to determine the frequency 
and extent of maintenance and overhauls without jeopardizing the 
integrity of equipment in satisfactory condition by disassembling it 
for visual inspection. In addition it warns of an impeding failure 
that could be hazardous to personnel and adjacent equipment. 

Vibration analysis has proved to be a powerful tool for analysis of 
failure in rotating machinery. Vibration analysis with appropriate 
instrumentation can provide reliable knowledge of machine condition 
and allows relatively inexpensive repair on a prescheduled controlled 
basis in comparison to catastrophic failures whose repair may cost 
many more times. Condition monitoring based on vibration signature 
analysis. is a custom built system for effective product quality 
assurance or machinery fault diagnosis. It normally includes the 
fol lowing- 
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1. Listing of all critical locations of a machine for monitoring the 
signatures. 

B. Identifying the vibration generating mechanisms and establishing 
acceptance level of machinery vibration and its spectral 
characteristics. 

3. Determining the machine condition and normal vibration level and 
spectrum . 

4- Selecting the interval for periodic vibration checks. 

5. Starting from a sample data recording system which correlates 

t 

the vibration signatures with machine condition. 

Advances in data acquisition techniques and the developments in the 
areas of knowledge achieving, processing and programming offer 
potential for application of automated strategies for machinery 
process supervision which would permit early detection of faults than 
is possible through conventional limits and trend checks. 

Meaningful and more efficient strategies for condition monitoring of 
rotating machinery are a pressing need , due to increasing demands on 
reliability and safety, especially for power generating systems. 

1-1 Literature survey 

Diagnostic systems generally use a multitude of observations for 
diagnosis and use diagnostic models to define the relationship of the 
measured signals to individual faults (Pau, 1975). Diagnostic models 
are either mechanistic or empirical. Mechanistic models which are 
derived from the physics of the process are dynamic state models and 
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assume complete knowledge of the measurement characteristics 
(Willsky, 1776 ; Ishermann, 1984). These models have an analytic 
foundation, therefore they can take advantage of the various 
parametric identification techniques for diagnosis (Stein and Park 
1988 ; Ono et al , 1987). 

Lipsom has presented a general discussion on the failure of machine 
parts. Elaborate results of experimentation on unbalance producing 
phenomena and balancing by statistical and probabilistic methods have 
been presented. Rao (1991) has presented a documentation of the 
dynamic performance of a machine due to rotating unbalance, fluid 
film bearings , asymmetric cross sections of the shaft etc. Most 
mechanical systems are, however too complex for such mechanistic 
modeling and the approach has limited application. 

Empirical models on the other hand are developed from heurist ical ly 
defined and / or experimentally established fault signatures. Since 
these models do not require detailed understanding of the process, 
they gain wider acceptance for performance monitoring of complex 
systems (Collacott, 1976 ; Compel, 1988 ; Lyon and Dejong, 1984 ; 

Gilmore and Gingher, 1986). 

Ignezio (1969) , Alty (1989) , Rich (1990) present the expert system 
building tools in their works. Hill and Beines (1984) painted out the 
utility of an expert system for maintenance personnel and suggested 
the development of one using Sahre‘5 database <1968). An expert 
system for machine monitoring is described by Reilly and MacPhail 
(1991). Lewis and Seth (1990) also described a rule based expert 
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system for vibration oriented diagnosis of turbomachinery for fault 
identification and predict maintenance in PC based Prolog 
environment- Typically the experts knowledge is represented by if — 
then rules and the observations are contained in a fact— base- An 
inference engine selects the rules to apply by matching the 
conditional side of the rules to the observat ions - 

l.E Present Study 

The present study concerns itself with development of diagnostic 
shell for rotating machinery with special reference to power plants- 
In the context of power plants very significant experimental 
investigations have been carried out by Sohre (1968)- A vast amount 
of data, for various kinds of faults in the rotating machinery and 
the associated symptoms in vibration signatures have been generated 
by Sohre for a real life power plant. An attempt has been made in the 
present study to organize Sohre *s data in the form of a knowledge 
base- Automation of the decision making process by relating 
symptoms to the cause through development of inference engines is 
carried out- Two separate inferencing techniques have been developed 
- one based on probabilistic reasoning based on Baye's Rule and the 
other based on Multi - Valued Influence Matrix ( MVIM ) approach 
suggested by Danai and Chin (1991), which in preference to 
probabilistic approach tries to match the symptom vectors from the 
test to those related to a specific fault, thereby achieving a 
reduction in decision making time-lag in comparison to a 


probab i 1 ist ic approach - 



A brief review of the commonly occurring faults in rotating machinery 
is presented in Chapter B while Chapter 3 reviews expert systems 
programming procedures- Chapter h presents the modeling of the 
problem and shell structure- Chapter 5 presents a case study on a 
power plant evaluation and Chapter 6 concludes with a discussion on 
scope for future work . 
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CHAPTER 2 


PROBLEMS IN ROTATING MACHINERY 

Rotating machinery are employed in various industrial installations 
such as steam and gas turbines , turbogenerators , internal 
combustion engines , reciprocating and centri'fugal compressors etc. 
The complexity of the dynamics of these machinery depends on the size 
as well as the operating speed. On account oT the ever increasing 
demand Tor power and high speed transportation , the rotating 
machinery is being made increasingly light and flexible, thereby 
increasing the dynamics complexities. Any fault in the system results 
in corresponding changes in the behavior of the system and vibration 
signatures of the machine become important indicators of the health 
of the system Typical rotating machinery installations involve 
mechanical elements like rigid or flexible shafts carrying one or 
more discs with bladed components. The arrangement may be supported 
in rolling element or fluid film bearings with associated components 
like couplings, seals, gear-boxes, pulleys etc. The installation is 
carried by pedestals which are mounted on foundations,. The commonly 
occurring faults in rotating machinery are cracks on shafts, loss of 
blades on discs, deteriorating balance, bows of shafting, shaft 
misalignment, oil film whirls, gear misalignment and damage, bearing 
eccentricities, coupling inaccuracies and damage, pedestal and 
foundation cracks etc. 
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e.l Classification of Rotating Machinery Vibration Causes 

The various kinds of rotating machinery faults become the cause for 
changes in the dynamics of the system. The changes in the system 
dynamics are reflected in the vibration signatures of the machine. 
The causes of vibration in the rotating machines are broadly 
classified as follows : 

(1) Unbalance 

Unbalances in machinery arise due to eccentricities 
between geometrical and mass centers in components. During 
rotation these unbalances generate large centrifugal forces in 
the machine which cause vibrations at the rotor speed i.e. one 
per revolution. Some specific unbalances created during the 
running of the machine give rise to additional characteristics, 
also Fig 2.1 lists causes of unbalance and the frequency of 
vibration they cause. 

(2> Distortion ( Casing and Foundation ) 

Casing and Foundation distortion .cause vibration in an indirect 
way either by generating misalignment between driver and driven 
machine or by causing internal rub or uneven bearing contact - 
This in turn transmits forces to the rotor, inducing it to 
generate forces of its own such as unbalance and a wide variety 
of oil film and friction induced forces. Another possibility is 
that the loadd on casing supports shift and can set off a 
series of resonance problems. Piping forces and foundation 
distortion often cause this type of difficulty. 
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Cause of unbalance 

Observable sign 

Frequency of 

vibration 

Disk or component 

eccentric an shaft 

De tec tab 1 e runou t on 

slow rotation 

One / rev 

Dimensional inaccuracies 

Measurable lack of 

symmetry 

One / rev 

Eccentric machining or 

forming inaccuracies 

Detectable runout 

One / rev 

Obi ique angled 

component 

Detectable angular runout 

One / rev 

Bent shaft . Distorted 

assembly 

De tec table r unou t on 

slaw rotat ion ; heavy 

vibration during rotation 

One / rev 

Section of blade or 

vane broken off 

Observable bearing 

vibration dur ing operation 

One f rev 

Eccentric accumulation of 

process dirt on blade 

Bear ing vibration 

One / rev 

Dif ferent ial thermal 

expansion 

i Shaft bends & throws c.g. 

out. Source of heavy 1 

vibration 

One / rev 

Nonhomogeneous component 

struc ture , subsurface 

voids in casting 

Rotor machined concent r ic . 

Bearing vibration during 

operation 

One / rev 

Non-uniform process 

erosion 

Bearing vibration 

One / rev 

Loose bolt or component 

slip 

Vibration reappears 

after balancing 

One / rev 

Trapped fluid inside 

rotor 

Vibration reappears 

after balancing 

One / rev 

Ball-bearing wear 

Bearing vibration 

One, two or 

higher / rev 


Table P.l 


Typical causes of rotor unbalance 

































(3) hi sal ignment 

Coupling misalignment causes friction or deflection forces in 
the coupling and this in turn causes the rotor and bearing 
system to deflect and to amplify these forces possibly creating 
secondary phenomenon such as harmonic resonance which can become 
severe^ Frequencies are characteristically E x Running 
Frequency (RF) for gear couplings, but 1 k RF is also often 
observed - 

(4) Bending Criticals and Resonance 

A rotor runs into resonance when the rotational speed is equal 
to its lateral (bending) vibration frequency* Bending critical 
speeds can be easily detected by synchronous whirl conditions 
and fairly large amplitude at the rotor speed* The critical 
speed region occurs over a fairly large range in the case of 
oil film bearings and is often accompanied by backward whirl. 

Resonance of the structures, support and auxiliaries cause 

/ 

•fairly large amplitudes of vibration at the rotor speed. 

(5) Asymmetric shaft 

The response of the asymmetric sha-ft has several harmonics and 
the frequencies observed are 1 X rev, E X rev, 3 X rev and 
sometimes even harmonics, if the symmetry is predominant. 

(6) Torsionals 

Torsional criticals in machinery are detected by torsional 
vibration measuring instruments. Unless there is a gear box in 
the transmission unit, it is difficult to detect these criticals 
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by vibrometers or accelerometers used for conventional 

non-angular vibratory measurements Torsional criticals are 

predominant in reciprocating machinery, as there is large 
excitation torque . In rotating machinery , they are not easily 
excited , unless there is a disturbance arising out of gear 
transmission units. Reciprocating machines have several 

criticals, some of them major, between the starting and 
operating speeds. In pure turbomachinery torsional excitation is 
normally very weak and inconsequential. 

(7) Oil Film Whirl 

Instabilities in rotors can also arise due to oil film forces in 
the bearings. This is self excited type of motion and is called 
oil whirl or oil whip. The frequency of response at the onset of 
the whirl is always a little less than 1/E rev , in the region 
of 0.4 to 0.5 X rev. 

(8) Mechanical Looseness 

Looseness of components is the likely cause of most number of 
problems. Loose rotor components such as disks, sleeves, thrust 
collars etc-' cause internal friction problems. Internal friction 
is a destabilising factor. The frequency of vibration is always 
the rotor critical speed and hence easy to detect. Mechanical 
loose components like bolts, give rise to 1 X rev and harmonic 
frequency signals due to secondary phenomena. The amplitude and 
phase continually change. Loose assembly of bearings give rise to 
subharmonic response and the typical frequency response is 
1/E X rev and 1/3 X rev. 


10 



Table E.E Bummarizes some of the common faults in the rotating 
machinery along with the resultant frequencies and directions of 
dom inant vibrati ons . 

8-E Comments 

The large variety of components on rotating machinery installations, 
with their complex dynamic interrelat ionsh ips and the possibilities 
of occurrence of different kinds of faults, make health monitoring of 
such a difficult task* Conventional manual trend checking techniques 
are obviously inadequate which has led to increased activity in 
development of automated diagnosis and prognosis procedures for 
rotating machinery installations- 
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Paul t 

1 Rota t i ng unbalance 
E Reciprocating unbalance 

3 Bent shaft or coupl ing 
^ Mi sal ignment 

5 Mechanical looseness 

6 Loose Rotor components 

7 Loose journal bearing 

8 Oil film whirl 

9 Damaged hydrdynamic 
bearing 

10 Assymetric shaft 

Abbreviations 

R Radial 

A Axial 

B Bearing housing 


Dominant Frequency Direct 

Rot at i ng speed , n R 

n , En , 3n R 

n/E y n/3 etc. 

n R , B 

En , sometimes R ^ A 

n y 3n , 4n B 

n , En , 3n R , A 

B 

Rotor critical speed R 

and higher harmonics B 

n/E , n/3 R, B 

0.4 to 0.5 R 

and harmonics B 

Change in higher R 

•frequencies B 

n , En , 3n R, B 


Table E-E 


Vibration Characteristics of Some Machine Parts 



CHAPTER 3 


EXPERT SYSTEM PROGRAMMING 

Expert Systems are 'intelligent' computer programs that exhibit 
within a speci'fic domain a degree of expertise in problem solving 
that is comparable to a human expert. Here intelligent refers to 
Artificial intelligence (AI). Expert systems are considered to be one 
of the most successful application of Al . Based on its function an 
expert system can be classified as being of the following type z 
1. Prediction System : Infers likely consequences from given 

situation .(eg. Expert system for weather-forecasting) 

S. Diagnostic System : Infers system malfunction from observed 

data / phenomenon (eg, the well-known Expert system for 
diagnosis of blood infection : MYCIN). 

3, Design system i Solves or assists to solve design 

problems and develop configuration of objects to meet specific 
design constraints of the problem (eg. PROCON for providing 
assistance for design of Mechanical systems. ) . 

Planning Systems : Designs a sequence of actions based on 

requirements and criteria (eg. Expert system for process 
planning in Arc Welding)^ 

5. Speech Recognition : Tries to interpret spoken input <eg. 

HEARSAY I and ID. 

6. Monitoring System : Compares the system performance with 

model of expected performance (eg Process Intelligent Control , 
PICDN for monitoring a plant operation control system ). 
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The motivation for building enpert system consists of one or more 
of the following : 

1. to act as substitute of the costly human expertise 
to formalize and standardize expert knowledge 
3. to integrate diverse source of knowledge about a problem domain 
and its solut ion methods - 

3.1 Organs of Expert Systems ' 

Some important terms associated with knowledge based systems are as 
f ol lows : 

* Knowledge base . — The knowledge base is the most important part 

of an expert system . It is in the form of facts and rules. The 
facts are termed deep knowledge. The rules are simply the 
heuristics which are also termed shallow knowledge- 

* Inference Engine - The inference engine serves as the 

inference and the control mechanism for the expert system and as 
such is an essential part of the expert system- It interacts 
with the user and uses the facts in the knowledge base according 
to the rules contained it ( Fig 3.1 ). Like the knowledge 
base the inference engine also contains rules and facts. But 
while the rules and facts of the knowledge-base pertain to the 
specific domain of expertise, the rules and facts of the 
inference engine concern more with the general control and 
search strategy employed by the expert system- 




User Knowledge 
engineer 


Fig, 3.1 A generic expert system 






ExpFrt 5y5te»m Shells ^ It is a commerc ial ly available expert 
system architecture and it has the main components of the 
expert system except for the knowledge base* The knowledge^-base 
is inserted into it to make a complete expert system. The use o-f 
shells thus frees the knowledge engineer from the need to 
repeatedly develop all the supporting elements of an expert 
system and thus to focus his or her attention on the development 
of the knowledge-base , Some well known commercial shells are 
Rulemaster, CxSYS, EXPERT etc. 

Working Memory - The contents of the working memory consist 
of facts. However unlike the facts within the knowledge base 
these facts are those that have been determined for the specific 
problem under consideration during the consultation session. 
More specifically the results of the inference process are new 
facts and these facts are stored in the working memory. 

Rule Adjuster - This serves as a rule editor, i«e«, it enters 
the rules into the knowledge- base during the development phase 
of an expert system- It mainly concerns with expert systems 
built on expert system shells. 

Knowledge Representat ion - The knowledge that is contained 
within an expert system consists of the following - 
A priori Arfo^tiogfge - the facts and rules that are known about a 
specific domain prior to any consultation session with the 
expert system. 

Inforretf AnaffioIgB - the facts and rules concerning a specific 



case that are derived during, and at the conclusion of a 
consultation with the expert system. Rules are generated during 
a consultation with an expert system only if the expert system 
is capable of learning. In general inferred knowledge consists 
simply of new facts, or conclusions. 

The major concerns of knowledge representation are the following 
To provide a format compatible with the computer. 

To maintain as close as possible a correspondence between this 
format and the actual rules and facts. 

To establish a representation that can be easily addressed, 
retrieved, modified and updated. 

Borne common modes of knowledge-representation are as follows 
Dbject-at tribute-value Triolets (OAV Triplets) 

OAV triplets provide a particularly convenient way in which to 
represent certain facts within a knowledge base and may be 
extended to provide the basis for the representation of 
heuristic rules. Using airplane as an example some of the 
attributes include the following : 
o Number of engines 
o Type of engine 
o Type of wing design 

Values may be symbolic or numeric. And in the above example 
the values may be as given below 
o Number of engines = 4 

o Type of engine = 
o Wing design = sf/ept bscA 
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Semantic Networks 


A semantic network may be thought of as a network that is 
composed of multiple OAV Triplets in network form. However 
rather than pertaining to just one attribute for a single 
object, semantic network can be used to represent several 
objects and several attributes per object. 

Frames 

This is a particularly robust way of representing knowledge. A 
frame contains an abject plus slots Tor any and all information 
related to the object. The contents of such slots are attributes 
and attribute values. However in addition to staring values for 
each attribute, slots may contain default values, pointers to 
other frames, and sets of rules and procedures that maybe 
implemented. Thus semantic networks may be thought of as a 
subset of the concept of frames. Frames are useful for the 
design of large““5cale, complex expert systems - particularly 
those involving a large amount of a priori facts and multiple 
ob jects . 

Representation via Logic Statements 

The most common form of logic is propositional logic. A 
proposition in turn is a boolean statement , it may either be 
true or false. Propositions may either be linked together with 
various operators called logical connectives such as AND, OR, 
NOT and EQUIVALENT. Linked propositions are called compound 
statements . 
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Predicate Calculus is simply an extension of propositional 
logic* The fundamentals of predicate calculus are the object and 
the predicate* A predicate is simply a statement about the 
object or a relationship that the abject possesses. Some 
examples of the use of predicate calculus are as follows : 
o Tail(dog) , which is read as 'a dog has a tail* 

o Brother(Ram, Shyam) , which is read as 'Ram is Shyam's brother’. 

Representation via Rule-Based Systems 

This is the most popular mode of knowledge representation ai: 
least at this time. The rules are of if-then type or production 
rules . 

Some of the important properties of these type of rules are as 
follows : 

o Name ( name of the rule > 
o Premise ( the if portion of the rule ) 

D Intermediate conclusion < the then portion of the rule ) 
o Conclusion ( the then portion of the rule ) 
o Notes (notes associated with the rule ) 

a Confidence factors ( a measure of confidence in the rule s 
conclusion 

o Chaining Preference ( the normal or default mode of search by 
rule 


th 
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3.E? Comments 


Few expert systems involve all these organs together. The type of 
knowledge base and inference engine depends on the problem 
concerned. The present expert shell uses lists of lists to represent 
the matrix of probabilities because the multiple array data structure 
is not available in GCLISP.lt comprises of a Knowledge Base, a Logic 
Analyzer, a Network and Inference Engines. 
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CHAPTER 4 


DIAGNOSTIC SHELL 

The diagnostic shell developed during the present study is based on 
empirical data by Sohre (196B) and available heuristic knowledge- As 
stated earlier, an empirical model does not necessitate detailed 
analytic understanding oT the process and thus finds greater 

acceptance in terms of efficiency. Since the knowledge or data in 
most cases may be incomplete or uncertain, empirical models employ 
probabilistic reasoning techniques such as Baye*s Rule, Fuzzy Logic, 
Dempster— Shafer Calculus etc. The present uses shell Baye * s rule for 
probabi 1 ist ic reasoning. In addition an Influence Matrix Approach, 
Danai and Chin (1991) for reasoning is also presented in this study. 

^.1 Shell Structure 

The diagnostic shell developed is described in Fig 4.1. The shell 
comprises of a Knowledge Base which is organized in the form of 
frames. A Logic Analyzer is appended to each frame For 

probabi 1 isation reasoning within the frame. The Logic Analyzers of 
the various frames are networked for overall probabilistic reasoning. 
User interface is available on initiating the shell whereby the 

observed symptoms can be logged in. At the network level user 

interface is once again provided to choose to an appropriate 
inference engine for d i agnos i s - 

The various components of the shell developed are described in the 
fol lowing^ 

central LiBRARX 

3, I I T.. KANPUR 


t-a *^1 



To Select 
Appropriate 
Inferencing 
Technique 



User 

Interface 




Machine 

History 

Frame 


Frame 

Logic 

Analysers 





Fig. 4,1 Shell Structure 




Knowledge Frames 


The knowledge base is organized in the form o'f “'frames*. A 
frame, in the shell developed in this endeavor, is a collection of 
probabilistic data relating 'causes* and 'symptoms* for a particular 
'monitoring parameter*. For example the frame in Fig 4-2 would 
contain the probabilistic data on the occurrence of particular 
'symptoms* namely vibrations at Rotor Frequency , Oil Whirl 
Frequency, 1 X Running Frequency, E X Running Frequency due to 
'causes* of Unbalance , Casing Distortion, Bearing Eccentricities. 
The symptom in the Fig 4.H are all related to the monitoring 
parameter of 'Frequency of Predominant Vibration* . 


CAUSiCSi or 
VIBRATION^ 

SYMPTOM 

ROTOR OH 
STATOR 
RESONANT 
FREO 

40-505^ 
of O IL 
WHIRL 
FREO. 

SO-- 

i K 

RUNNINO 

FREO. 

2 X 

R. r. 

H 1 OHER 
MULTI- 
PLE 


0. 05 

O. C 

O. O 

rt jjj 

O 0 5 

0 . 05 

CASINO DISTORT ION 

D 

O 

10 

C. 04 

0.03 

c. <s 

O 2 

4 

BRG-ECCENTRICI TY 

O. 

O. O 

O. O 

0. 4 

O. O 

1 o o 


Fig 4,S Cause - symptom frame for monitoring parameter 

' Predominant Frequency * 


Fig 4-3 shows a frame for the same set of 'causes* but 'symptoms* for 
monitoring parameter 'Direction of Predominant Amplitude*. 


CAUSES or 
VIBRATION^ 

SYMPTOM 

VERTICAL 

HOHIE. 

AXIAL 

UNBALANCE 

O. 4 

O. 5 

O . 1 

CAS INC DISTORTION 

O. 4 

O. 5 

O . 1 

BRC -ECCENTRIC I TY 

O. 4 

O . S 

a . 1 


Fig 4.3 Cause - symptom frame for monitoring parameter 
' Location of Predominant Amp 1 i tude * 

\ 


























In the shell developed the ‘frames' incorporate the 
'causes' (causes are the rows in frames) s 

( 1 ) Unbalance 
(E) Permanent Bow 
(3) Temporary Bow 
{^) Casting Distortion 

(5) Foundation Distortion 

(6) Leaf Cut 

(7) Rotor Rub 

(8) Misalignment 

(9) Piping Losses 

(10) Bearing Eccentricity 

(11) Bearing Damage 
(IE) Excited Vibration 
(13) Unequal Stiffnesses 
(1^) Thrust Bearing Damage 

(15) Loose Assembly of rotor 

(16) Gear Inaccuracy 

(17) Coupling Inaccuracy 

(18) Aerodynamic Excitation 

(19) Rotor System Critical 
(SO) Coup 1 ing— critical 
(El) Overhang-critical 

(EE) Structural Resonance of Casting 
(E3) Pressure Pulsations 
(S4) Electrically Excited Vibrations 
(E5) Vibration Transmission 
(E6) Oil Seal Induced Vibration 


E-ff 
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The frame incorporates seven frames in its knowledge base. Each frame 
pertains to a specific ‘monitoring parameter*. The seven monitoring 
parameters along with the associated symptoms (which form the columns 
in the frames) are : 

( 1 ) Predominant Frequencies 

Rotor or stator resonant frequency 

40 - 50 '/. oil whirl frequency 

50 - lOO’/i oil whirl frequency 

Equal to running frequency 

Twice running frequency 

Higher multiple of running frequency 

Half of running frequency 

Quarter of running frequency 

Lower multiple of running frequency 

An odd frequency 

Very high frequency 

(H) Direction of Predominant Amplitude 
Vertical 
Hor i Eontal 
Axial 

(3) Locat ion of Predominant Ampl i tude 
Shaft 
Bearings 

\ 

Casting 

Foundation 


E5 



Piping 
Coupl ing 

(4) Amplitude Response Speed Variation 
^ Speed Com i nq Up ) 

Stays same 
Increases 
Decreases 
Peaks 

Comes suddenly 
Drops out suddenly 

Ampl i tude Response to Speed Variation 
^ Speed Slowing Down ) 

Stays same 
' Increases 
Decreases 
Comes suddenly 
Drops out suddenly 

(6) Predominant Sound During Vibration Test Run 
Low frequency rumble 
Hum 

Beat ( i e comes and goes periodically ) 

) 

High pitched whine 
Very high loud scream 
Very high squeal 

Ultrasonic (Measured with sound analyzer) 
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^ ) R esponses f or the effect of pperat i no cond i t ions 


Comes in or goes out at 
No effect 
Full load in 
Full load out 
No load in 
No load out 
Part load in 
Part load out 
Start up only 
After load dump or surge 

When other machine starts or shuts down 

Known operational evidence 
Seals rubbed 
Shaft bent 

Thrust bearing damaged 
Bearings wiped out 
Bearings fatigued 
Bearings babbit squeezed out 
Case distorted or cracked 
Misalignment 

Coupling burned or pitted 
Gear teeth broken or pitted 
Gear teeth marked on the back side 
Shaft cracked or broken 

Bailing or fretting marks under disks or hubs 
Coupling bolts loose 
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Foundation cracked 

Soleplate loose or rusted 

Sliding surfaces binding 

Thermal expansion restricted 

Fluid marks on the internals 

Rotor components eroded 

Solid accumulated on vanes and rotor 

Salt deposited on internal 

Flanges leaking 

Seals leaking 

The History of Machine 
Initial start-up 
Start— up during the first year 

Start-up after one year but before ten years 
Start-up after ten years 

A^-3 Probabilistic Logic Analyzer Using Baye’s Rule 



Fig Conditional Appearance Symptom 

Let Si. represent the occurrence of the ith symptom which may 
caused by any of the independent causes Ci, C 2 , Ca etc. 


be 
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If a symptom S is observed then the object of the Logic Analyzer is 
to predict the probability P < Clc / S ) of the existence of the 
cause Ck given the observation of the symptom 

With reference to Fig ^.3 the probability element pLj in ^ any 
knowledge frame can be seen as the percentage of cause Cj overlapped 
by symptom Si- 

Let S be the symptom of 'Foundation Vibration* in the 
frame concerning 'Location of Predominant Amplitude*, then this fact 
can also be represented as 

S = (Not Shaft) and (Not Bearing) and (Not casting) 

and (Foundation) and (Not piping) and (Not coupling) 
or 

B = (--, Si) f) (*n 52) p, (-n S2) p (S^) p (n S5 ) p (**-. Si5) (1) 

The knowledge frames consist of elements pLj = the probability of 
occurrence symptom Si when fault Ck is present in the machinery. 

If the appearance of symptoms is independent of each other, then 
P (Si p Sk) = P (Bl) * P (Sk) 

Application of O) in (1) gives 

P (S / Ck) = P (nSi / Ck) * P {HBz /Ck) * P(nBa /Ck) * P<S4 /Ck) 

# P (“S!S / Ck) * P (-nSiS / Ck) <3) 

Since P (nSL / Ck) = 1 - P (SL / Ck) <^> 

Using the available data in the knowledge frame P (S / Ck) can be 

evaluated. 

Defining the total probability, P(S), of occurrence of a symptom S 
as 
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P (S) 


k 

z p (B U Cl:) 
k 

2] P <S / Ck) * P ( Cl. ) (5> 

Where P < Cl. ) is the probability of the presence of a cause Cl;. The 
probability P ( Cl. ) is picked from the Machine History Knowledgt 
l-iame which is constructed as shown in Fig 4.5. 


CAUSES OF 
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FIRST 
YEAR . 
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A . O 

O. O 

O. O 

O. O 

PKRWANENT BOV 

O. 33 

o. o 

O. 33 

a. 3^ 

TEMFORARV BOW 

Cl. <S7 

Cl. 33 

O. o 

o, o 


Fig 4.5 Machine History Knowledge Frame 


Using the definition of conditional probability 
P(Ck/B) = P(Ckp, B) / P(B) 

Since 

P ( S / Ck ) = P ( B p Ck ) / P ( Ck ) 

or P { S n Ck ) = P ( S / Ck > ♦ P ( Ck ) 


( 6 ) 


(7) 


As P ( S fj Ck ) = P ( S p Ck ) , (6) gives 

p ( Ck / S ) == P ( S / Ck ) ♦ P < Ck ) / P ( 5 > <8) 

Employing (5) in (8) the conditional probability of the presence of 
cause Ck given the observation of symptom B can be written as 

P ( S / Ck ) * P ( Ck ) 


P ( Ck / B ) 


(9) 


J] P (B / Ck ) * P ( Ck) 


Equation (9) is known as Baye's Rule for Conditional Probability, 
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4 .^ 


Networking of Frames 


After obtaining the conditional probability of the existence of a 
cause given the observation of a particular symptom for each of 
monitoring parameter , the probability values for all the frames are 
to be networked to give an overall conditional probability of the 
existence of a cause for a given set of observations- The networking 
scheme is described in the following. 

The Odds* of the presence of a cause is defined as the ratio of the 

probab i 1 i t i es of occurrence of the cause to its nan— occurrence . For 

1 

example for an event X with the probability of occurrence as P (X) 
the odds are defined as 

0 (X) = P ( X) / C 1 - P (X) 3 (lO) 

Odds Likelihood Ratio of the existence of a cause Ck is defined as 

X = P ( S / Ck> / P ( S / -1 Ck) (11) 

Simplified the above gives 

X = □ ( Ck / B ) / 0 ( Ck ) (IE) 

The eight monitoring parameters are independent of each other. 
Therefore the overall likelihood ratio (X) is 

Xe = Xi # Xz ♦ Xa * X4 * Xs * Xd * Xv * Xb (13) 

This effective value Xe can be used to calculate the final 
probability of the cause 

Xe = 0 ( C / S ) / 0 ( C ) (1^5 

Where S is the global set of symptom observed pertaining to the 
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various knowledge frames. 

Therefore 

0 < C',. / S ) = 0 ( Cl. ) 

From this P ( Cl. / 5 ) can be calculated as follows 

P(Cl./B) = 0(C1. / S) / Cl+0(Cl. / S)3 

Thus P ( Cl. / S )can be calculated for each of the defects 
C4 etc and from these the defects with higher probabilities 
short listed. 

4.5 MVIM Method 

An alternative inference engine based on Multi Valued Influence 
Matrix diagnostic reasoning is also built into the shell. The MVIM 
method involves matching the vector of observed symptoms against 
vector of symptoms pertaining to a particular cause. The inference 
based on MVIM reasoning is more suitable in cases when data in a 
particular knowledge frame is uncertain or doubtful and no data 
entries are made in the frame grid. The absence of the would be taken 
zero probability leading to erroneous ordering of the causes. The 
MVIM method works across the frames and ignore the absence of data. 

As stated diagnostic reasoning in the Multi Valued Influence Matrix 
method is based on matching the vector of observed symptoms against 
vector of symptoms pertaining to a particular cause The particular 
influence vector at the closest angle to the vector of normalized 
measurement vector is considered to correspond to the occurred fault. 
Thus every fault Ck may be characterized by its influence vector as 


( 15 ) 

< 16 ) 

Cz, Ca, 
can be 
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follows 



where at = Prob (Sv / Ck ) 

< Probability that symptom Si is observed 
given fault Ck ) 
thus O ^ < 1 

Let the observed vector be 



It is to be noted that Y is a binary vector, so that vi taking a 
value O or 1 , 1 if ith symptom is observed and O if it is not 
observed . 


Since the angle & between two vectors 


A = 

Aii 

+ Azj 

+ Ask 




B = 

Bii 

+ Bzj 

+ Bsk 




is 







6f = Inverse 

cos 

(Ai Bi + 

Az Bz + 

Aa B3> / 

A . B 

{ 17) 

The angle between 

the 

observed 

V vectors 

and the 

symptom vector 

D 

> 


cause Ck can be found out by extending eq.(17) to n dimensions. 

As an example consider the symptom vector for the cause unbalance 
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Where 0.6 is the probability that vibration occurs at the 
shaft . 

Let the observed vector of symptom be 

■ 0 
1 

V = 0 

0 
0 
0 

Therefore 

Vk = y ( 0.9 » 0.9 + 0.1 * 0.1 +0+0+0+0) 

= y O.BE 
= 0.905 
Y = 1 

y = Y 

thus C.OB& = 0.1 * 1 / 0.905 

= 0.1105 

The cosine of the angle between the each of the cause vector can thus 
be calculated and the causes which are at lower angular separation 
to the symptom vector can be short listed. 
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i^.b Comments 


The shell is developed in GCLIBP. LISP stands Tor LISt Processor and 
GCLISP is abbreviation for Bolden Common LISP marketed by Golden 
Hill Corporation and is a small memory PC based version of LISPUnlike 
the procedural languages (like Fortran, Cobol , C and Basic) ; Lisp 
and Prolog are object oriented languages and they are more suited 
Tor developing expert systems and AI. While Prolog employs object 
data and their relationships, emphasizing symbolic processing, LISP 
is more suited for backward chaining operations. 
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CHAPTER 5 


SHELL USAGE 

The usage of the Expert shell is described in this chapter. Diagnosis 
of an high speed turbomachinery installation is illustrated. The 
knowledge base is constructed on the basis of the exhaustive data 
available from Sohre's (196B) work on Petroleum Plants. Employing the 
various features of LISP programming a question-answer based 
interface is developed. The output from the shell is obtained in the 
form of most prominant faults likely to be present in the machinery. 

5.1 Knowledge Base Creation 

Typical operational knowledge for a High Speed Turbomachinery 

installation in a Petroleum Plant is organised in the form of 
‘Frames'. Each frame pertains to a specific monitoring parameter. 
Table 5.1 to 5.0 describe the organisation of these Cause-Symptom 
frames for various monitoring parameters . The grid values represent 
the Cause-Symptom relationships in the form of probability of 
occurence- The probabilistic knowledge is fed to the computer in the 
form of list of lists, for example the symptom Location of 
Predominant Amplitude the list of list would of the following form 


( 0.9 

0. 1 

o 

b 

o 

o 

0.0 

0 . 0 ) 

( 0.9 

0 . 1 

0.0 

0.0 

0.0 

0 . 0 ) 

( 0.9 

0 . 1 

0.0 

0.0 

0.0 

0 . 0 ) 

( 0.9 

0. 1 

0.0 

0.0 

0.0 

0 . 0 ) 
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Table 5.1 s Machine History Knowledge Frame 
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Tablp,5.3-: Monitoring Parameter Direction of 
Predominant Amplitude 
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Table 5.5 : Amplitude Response to Speed Variation During Vibration-Test F 
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Table 5.7 s Manitaring Parameter Predominanat Sound (During Vibration Test Run) 


















































































































































































































fable 5-8ct : Manitaring Parameter Known Operational Evidence 























































































































































































































Table 5.8b : rianitorinq Parameter Known Operational Evid 



































































































•Car' and 'cdr* operations are done to feed / alter elements of the 


list. 

5.S User Inerface 

User interface is provided to feed in the observation on the machine 
under diagnosis. On execution the shell progresses frame by frame and 
puts forth a set of questions to the user for his / her responses. A 
typical question-answer session with the user is described below. 

(The output of the program is given in 'standard font' below while 
the user responses are in the italics font- Explanation is starred 

(»> . ) 


WELCOME TO THE PROGRAM FDR DETERMINATION 
OF FAULTS IN ROTATING MACHINERY 

SCf{^£N / 


SELECT THE LOGIC ANALYZER 

Do you want to use the BAYE's method (Type 1) 

Da you want to use the MVIM method (Type E) 

Type any other key to exit 
/ 


SCf(££fi £ 


A6 





I 

Please answer the following in yes / no only . 

If you chance to feed in incorrectly somewhere then carry 
on, you can correct in the end. 


S£/i££(V 3 


FRAME FDR PREDOMINANT FREQUENCIES 

Is the predominant frequency equal to rotor or stator resonant 
frequency 

Is it equal to 40 - 50 Vi oil whirl frequency 
Is it equal to 50 - lOO X oil whirl frequency 
Is it equal to running frequency T.-.fes 
Is it equal twice running frequency 1 . . .fl/a 

Is it a higher multiple of running frequency ? 

Is it half of running frequency T.-.lVc 
Is it quarter of running frequency 

Is it a lower multiple of running frequency T.-.Aiff 

Is it an odd frequency 

Is it a very high frequency 7.,. No 


FRAME FOR THE DIRECTION OF PREDOMINANT AMPLITUDE 

Is 

the direction 

vertical 7 .. .No 

Is 

the direction 

horizontal ?. ■ mVes 

Is 

the direction 

axial ?-- - 


SC£££N 5 
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FRAMF 

1 

FOR 

THE LOCATION OF PREDOMINANT AMPLITUDE 

\ 

Is 

it 

in 

the 

shaft ? . . .Tes 

Is 

it 

in 

the 

bearings ? * . ,^^<7 

Is 

i t 

in 

the 

cast ing ? - . .A^£7 

Is 

it 

in 

the 

f oundat ion ? . , - 

Is 

i t 

in 

the 

piping ? . , .da 

Is 

it 

in 

the 

coupl ing ? . . .fl/o 


S£;/l££N 6 





FRAME FOR AMPLITUDE RESPONSE 10 SPEED VARIATION 
< SLO WING DOWN ) 

Does it stay same ?...(Vc' 

Does it increase 

Does it decrease l...yes 

Does it come suddenly 1 . . .fllo 

Does it drop out suddenly 

SC/{££IV B 


FRAME FOR PREDOMINANT SOUND 
< DURING VIBRATION TEST RUN ) 

Is the predominant sound a low frequency rumble 
Is it a loud roar ?..-A'es 
Is it a hum ?...No 

Is it a beat ? (i e Does it come and go periodically ?) 

Is it a high pitched whine ? ffo 

Is it a very high loud scream 1 . . .No 

Is it a very high squeal 1 . . .No 

Is it ultrasonic (measured with sound analyzer) 7... No 


. .No 


S££££N P 




FOR RF^P^SEB FDR THE EFFECT OF OPERATING CONDITIONS 
COMEG IN OR GOES OUT AT 

No effect l...yps 
Full load in 
Full load out 1 . . .No 
No load in ? . . .A'/? 

No load out ? . « .fifo 

Part load in 

Part load out ? . , .No 

Start up only , .No 

I 

After load dump or surge ? , . .No 

When other machine starts or shuts down 1 . . .No 

SCREEN 10 


FRAME FOR KNOWN OPERATIONAL EVIDENCE 

Are the seals rubbed ‘I...Ybs 

Is the shaft bent I...N 0 

Is the thrust bearing damaged 1 . . .No 

Have the bearings wiped out 1 . . .No 

Have the bearings fatigued '{...No 

Have the bearings babbit squeezed out "{...No 

Has the case distorted or cracked ? No 

Is there misalignment ? No 

Is the coupling burned or pitted ? No 

SCREEN // 
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FRANr rOR KNOIJN OPERATIONAL EVIDENCE 

Are the gear teeth broken or pitted 

Are the gear teeth marked on the back side?... (Vo 

Ie> the shaft cracked or broken 7... No 

Are there gailing or fretting marks under disks or hubs ?.../Vo 

Are the coupling bolts loose I... No 

Is the foundation cracked 1 . . .No 

Is the soleplate loose or rusted 1 . . .No 

Are the sliding surfaces binding 1 . . .No 

Is thermal expansion restricted 1 . . .No 

Are there fluid marks on the internals 1 . . .No 

Are the rotor components eroded ? . . .No 

Has solid accumulated on vanes and rotor ?...f£’s 

Has salt deposited on internals l...yes 

Are main flanges leaking ? — No 

Are the seals leaking ? No 

S£/l££N // 


frame for the history of machine 

Has the machine had initial start-up ?...Xe5 
Has it had start-up during the first year 1 . . .No 

Has it had start-up after one year but before ten years I... No 
Has it had start-up after ten years 1 . . .No 

S£/I££N Ja 
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MODIFICftTTriN 


Da you want to faad again in the set of responses for any symptom 
(yes or no) , .IVa 

if ifSPr had printed yes fad Japing pauld have been printed on tfje 

screen 

Wr i te - 

1 “If you want to to refeed for predominant amplitude 
p „ Pqp direction of predominant amplitude 
3 “ for its Location 

A - For Amplitude Response - Speed Coming up 

f) “ For Amplitude Response - Speed Slowing down 

h _ Predominant Sound (During vibration test run) 

7 “ For Effect of Operating Conditions 

0 “ For Known Operational Evidence 

9 -- For History of Machine 


) 


SCHJBti 13 

5.3 Shell Output 

The shell output comprises of the list of the top seven defects in 
the order of decreasing probability. 
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DIAGNOSIS 


( Permanpnt~Bow 0.99911 > 
(Unbalance 0.99369) 
(Hearing-Damage 0,0.00183) 
( Temparary-Bow 0.0) 
(Casing-distortion 0.0) 
(Foundation-Distortion 0.0) 
(Seal— Rub 0.0) 




Do you want to know the MVIM method answer also. 

If yes type yes else type any other key to exit.... fas 

SCREEN JS 


DIAGNOSIS 

(Permanent-Bow 0.9055) 
(Unbalance 0.8756) 
(Temporary-Bow 0.733) 
(Bearing-Damage 0.7S08) 
(Foundation-Distortion 0.6A1A) 
(Overhang-Critical 0.H7E) 
(Leaf-cut 0.6S) 

SCEEE^ 16 
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g.4t Commpnts 


Thf shell is rtevelnpert in ^uch a „ay that additional fault and 
symptom-sot can bo added nithout much mod if icat ions . By changing the 
set or cau..,.., and th,. sot of symptoms the program can also be 
SKtendesd to any caso whore two parameters are related by cause-effect 

relationship . 



CHAPTER 6 


CONCLUSION 

A *,hf.n har, I, Pen developed and illustrated for diagnosis of 

fail] in in rotating machinery installations. Techniques of expert 
system programming have been employed for the development of the 
she ] 1 . 


The sfiel 1 (omprisc-s of a Knowledge Base in the form of Frames, a 
Prnbahi 1 i nt 11 . logir Analyzer and a Network. Two different engines, 
one bast'd on Ilaye's Method for probabilistic reasoning and the other 
based on Influence Matrix Approach are constructed for inferencing. 
The shell is illustrated for set Knowledge-Base for the case of a 
High Spf=>eri Turbomachine in a Petroleum Plant. 

Though tht* sliPll is developed with the view of specific application 
to rotating machinery, minor modifications may permit the usage of 
the shf=»ll to more general classes of machinery. The User Interface 
can be refined to give the user to define his/her own set of 
Monitor ing parameters and therefore the Knowledge Frames. The Causes 
and symptoms can be defined through the user-interface. Interface can 
be furiht'r can be further modified to allow for on-screen creation of 
knowledge base by the user. 

Feasibilities studies on the usage of the shell to cases like medical 
diagnosis and weather forecasting can also be undertaken. 
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